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We have analyzed the B° — > 7r°7r° puzzle in three kinds of models beyond the standard model 
(SM). It is shown that the minimal flavor violation (MFV) models, the minimal supersymmetric 
standard model (MSSM), and the two Higgs doublet models (2HDM) I and II can not give an 
explanation of the B° — > tv°tv 0 puzzle within la experimental bounds and the model III 2HDM can 
explain the puzzle without a conflict with other experimental measurements. If the constraint on 
Cs g from b —> s g is not imposed, for all kinds of insertions considered there are regions of parameter 
space, where the scalar quark mass is larger (much larger) than the gluino mass in the case of LR 
or RL (LL or RR), in which the puzzle can be resolved within lcr experimental bounds. 


The branching ratios of B° decays into two pions have been recently observed 

Br{B° -> 7 r° 7 r°) = (1.45 ± 0.29) x 10" 6 , 
Br{B° -> 7 T+ 7 T-) = (5.0 ± 0.4) x 10 -6 . 
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The large branching ratio of the B decay into neutral pion final states is unexpected. The decay amplitudes of 
B° —> 7 T 7 T can be generally parameterized as 
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where T,C,P, and Pew are the tree, color-suppressed tree, penguin, and electro-weak penguin (EWP) amplitudes 
respectively, and a = arg (^ 17 ) is the weak phase, where A p = V p bV* d (p = u,c,t). In SM one has the counting 
rules: the color-suppressed tree and penguin amplitudes are suppressed by a factor of A (A ~ 0.22 is the Wolfenstein 
parameter) and the EW penguin is suppressed by a factor of A 2 , with respect to the tree amplitude 3 . So one should 
expect by the naive counting rules that the branching ratio of the B decay into neutral pion final states is 0(A 2 ) 
of that for charged pion final states. However, the data (Eq. (1)) indict that the former is 0(A) of the later. The 
observed branching ratio of B° —> 7 r° 7 r° = (1.45 ± 0.29) x 10 ~ 6 is much larger than the theoretical prediction, about 
0.3 x 10 -6 , up to the a s order in the BBNS approach (QCDF) 33 in SM. In Li et al’s approach (PQCD) the leading 
order (LO) prediction ~ 10 ~ 7 3 the same order as that of the QCDF prediction. In the recent paper 3 the next 
leading order (NLO) PQCD calculations have been carried out and the results are that the ttK puzzle, the expected 
relation Hcp(H ± —> 7 t°A' ± ) « Acp(B° —> n ± K 11 ) disagreed significantly with the data, can be resolved but the 
predicted branching ratio of B° —> 7 r° 7 r° is about 0.3 x 10 ~ 6 which is still much smaller than the data, i.e., the 7 r° 7 r° 
puzzle remains. If the large branching ratio persists it could indicate new physics. 

Though B° —> 7 r° 7 r° is not a pure-penguin process and has the contributions from tree operators, the tree contri¬ 
butions are of the order same as the penguin contributions because of the almost completely cancellation between 
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the two terms in C 2 + C\/N c where 6% 2 are Wilson coefficients of tree operators, so B° —> 7r°7r° is sensitive to new 
physics. Therefore, it seems that a lot of new models beyond SM could enhance the branching ratio and consequently 


resolve the puzzle |8]. However, any new model must simultaneously give an explanation for the branching ratio of 
B° -> K°K° since the two processes are closely related at quark level: the flavor changing neutral current b —> d 
transition controls B° —> K°K° and the same transition gives significant contributions to B° —> 7r°7r° which are of 
the order same as the tree contributions in SM. Recently the branching ratio of B a —> K°K° has been measured as 
(0.96 ±0.25) x 10 
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|, which is consistent with the prediction from both the QCDF |4| and PQCD approaches 
Therefore, new physics (NP) contributions must satisfy the condition that they make Br of B° —> 7r°7r° enhanced 
but keep Br of B° —> K°K° basically unchanged, compared with those in SM respectively, which will impose the 
significant limit on NP models. 

In the paper we search for new models beyond the SM which can account for the data of branching ratios for both 
the B° —> n°TT° and B a —> K°K° processes. To be specific, we concentrate on the well-known three kinds of models: 
the minimal flavor violation (MFV) models, the two Higgs doublet models (2HDM) and the minimal supersymmetric 
standard model (MSSM). 


The effective Hamiltonian relevant for the two processes in the SM can be expressed as 
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where A p = V p bV* d , Qi ,2 and Qi(i = 3, ...10) are the tree and penguin operators respectively. Explicit forms for 
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C'i, C 2 , C'i, cv 7 , C$g and Q\, Q%, Qi, Q 77 , Q$ g can be found, e.g., in Ref. 

In the QCD factorization approach, the dominant contributions to the decay amplitudes are given by: 
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where the definitions of the parameters at and the chiral enhancement factors C, r K can be found in Ref. |3f- We take 
the values of running masses in the MS scheme for light quarks such that r" = r K = r hereafter. The electro-weak 
penguin and annihilation contributions are neglected in above formula, which leads to the 10% theoretical uncertainty. 

First we consider the MFV models. The MFV models beyond the SM discussed in the paper mean a class of models 
in which the general structure of flavor changing neutral current (FCNC) processes present in the SM is preserved. 
In particular, all flavor violating and CP-violating transitions are governed by the CKM matrix and the only relevant 


local operators are the ones that are relevant in the SM B- New parameters in the MFV models, e.g., the masses of 
charginos, squarks, Higgs particles in the MFV scenario of the MSSM, enter into Wilson coefficients of relevant local 
operators. Therefore, in the MFV models the amplitudes of the two decays are given same as Eq. Q (with values of 
a,i generally different from those in SM). 

We can model-independently determine \z\ = \ c [(°4 + raf)A p ] | from the measured branching ratios of B° —> 

7r°7r° and B° —> K°K° and a 2 = (0.0502 — 0.0689i)±(0.0025±0.0035*) which comes from the known tree contributions. 
The result is given in Fig. 1 where the 60% theoretical uncertainty (coming mainly from non perturbative parameters 
such as form factors, distribution amplitudes and CKM matrix elements) has been taken into account. There is a 
narrow region which can simultaneously fit the data. However, assuming Wilson coefficients of relevant operators, 
except the chromo-magnetic dipole operator, change a little compared with SM, which is the case in the MFV models, 
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FIG. 1: The constraints on 2 from the 
from B° —> 7r°7r° and the small lattice 



for the constraint from B° —> K°K°. 


K°K°. The big lattice denotes the constraint 


z in the region corresponds to 
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which can not be reached in the MFV models [10|,[11|. That is, the MFV models are excluded within lex experimental 
bounds. 

Next we consider models in which there are new operators in addition to those_in the SM, e.g., the 2HDM and 
MSSM. The effective Hamiltonian in the 2HDM and MSSM can be written as 
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where Q['\ i=ll,...,16, are the neutral Higgs penguin operators and their explicit forms can be found in Refs, 
with the substitution s —> d. The primed operators, the counterpart of the unprimed operators, are obtained by 
replacing the chirality in the corresponding unprimed operators with opposite ones. 

^From the effective Hamiltonian, Eq. ©, it follows that the main contributions to the decay amplitudes from the 
SM and new physics are given by: 
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where we have set rrid = 0. Due to the renormalization group equation (RGE) running, Wilson coefficients Ci, 
i=14,15,16, are related to C 13 and the known constants hi t 2 represent the running effects. The largest contributions 
to the hadronic elements of the neutral Higgs penguin operators at the a s order arise from penguin contractions with 
b quark in the loop, which are the same for B° —> and B° —> K°K° and have been included in 04 (see, for 

example, Ref. [0). Therefore, although they can enhance the branching ratio of B° —> 7 r° 7 r°, they alone can not 
resolve the puzzle because the branching ratio of B° —> K°K° will also be enhanced by them, which will not agree 
with the data. 

The new physics contribution, the terms proportional to C 11 and C 13 respectively, which contributes to one mode 
but not to the another (precisely speaking, the contribution to the another mode is TOd/m s suppressed) comes from 
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FIG. 2: The correlation between \Cs g (mw) — C' &g {mw)\ and \Ci 3 (mw) — C' 13 (mw)\- 


the hadronic matrix elements of Higgs penguin operators at the leading order in a s . It is the contribution that gives 
a possibility to make the data be account for without a conflict with all relevant experimental measurements. The 
key point is if one can have a sizable Cy} (/z) and/or Cy (/z) in the 2HDM and MSSM. 

Let us analyze how large C [3 and/or Cy are needed to fit the data. Let 

z= X ( a 4 + ra l) > 21 = ^02, z 2 =—[hi(Cu - C' n ) + h 2 {C 13 (n) - C' 13 (n))\, 
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From the data, 1.16 x 10 ~ 6 ^ Br( 7 r° 7 r°) ^ 1.74 x 10 -6 and 0.71 ^ Br(/?°/2°j ^ 1.21 x 10 _6 [3|i we have ri 2 > rn > 
r 22 > z' 21 . To satisfy the above two relations, it is necessary to have 


U2 + zi \ ^ rn - r 22 


( 10 ) 


In the model I and II 2HDMs and MSSM the Wilson coefficients of QCD penguin operators are not changed 
signific antly , compared with those in SM, and the Wilson coefficient of chromo-magnetic operator can have a significant 
change [16J. Taking the SM values of Wilson coefficients of relevant operators but the chromo-magnetic operator 
and using RQE running, we can obtain the correlation between \Cs g (mw) — G' Sg (m\y)\ and |Gi 3 (tow) — CJ 3 (wff)| 
from \z — Zi\ ^ rn and Eq. cm, which is shown in Fig. |2] where Gn = G 13 has been assumed for simplicity, 
without losing the generality of discussions*. It follows from the figure that \Cs g (mw) — C' 8g (mw)\min = 2.6 when 
Ci 3 (mw) — C[ 3 (mw ) = 0, which reduce to Eq. 0 in the MFV models, as it should be. 

It is well-known that the experimental upper bound of branching ratio for B s —> /i + ^z _ constrains severely pa¬ 
rameters in the MSSM and model I and II 2HDMs |n]] . Similarly, we show that the the corresponding bound for 


* In the figure \C% g {rnw) — C% g (rnw)\ has no upper bound because we do not impose \z — z±\ ^ r\ 2 - In all models beyond SM known 
so far, \C$, g (mw ) — C& g {m\v)\ never reach very large value (say, 5) when all relevant experimental constraints are imposed. We do not 
need to know the upper bound for the analysis in the paper. 
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Bd —-> fJ* implies that the Wilson coefficients of new operators in the MSSM and model I and II 2HDMs can not 

be large. The branching ratio Bd —■> fi + in the 2HDM and MSSM is given as 

Br(B d -> fj, + n~) = G 6 F 4 lT m k T B d fk\MWl- 4m 2 [(l - 4 m 2 )\C Ql (m b ) - C' Qi (m b )\ 2 + 

\C Q2 {m b ) - C' Q2 (m b )+ 2m(C w (m b ) - C' w {m b ))\ 2 ] (11) 

where to = In the moderate and large tan (3 cases the term proportional to (Cio — C' 10 ) in Eq. 111(1 can be 

neglected. The new CDF and DO combined experimental upper bound of Br (Bd —> g + g~) is 3.2 x 10 -8 jl^J at 90% 
confidence level. We have the constraint 


\J\C Ql {m w ) - Cq^mw )] 2 + I CQ 2 ( m w) - Cq 2 (tow )| 2 < 2.2 


( 12 ) 


where Cq| 2 are the Wilson coefficients of the operators Q['\ which are Higgs penguin induced in leptonic and 


semileptonic B decays and their definition can be found in Refs. [3Q- By substituting the quark-Higgs vertex for 
the lepton-Higgs vertex, it is straightforward to obtain Wilson coefficients relevant to hadronic B decays in the MSSM 
and model I and II 2HDMs. To translate Cq 1 2 into Cq x1 13 , we have Cq^ 13 {mw) ~ 0.037. Then it follows from 
Fig0 that \Cg,g{mw) — C 8 g {mw) | must be larger than 2.4 in order to resolve the puzzle. 

The Wilson coefficients C 8g in the b —> d transition are constrained by Br(H —> Xdg). Because there is no data 


X s g decays in the SM comes 
s g. In the presence of new 


for Br of the B —> Xd g decay and the difference between the B —> Xd g and B —* 
from CKM matrix elements, we assume the constraint on ci'l same as that from b 

“ * n 

physics a model-independent analysis gives that \C$ g (niw) — C$ (mw)\ < 2.01 when Br(b —> s g) < 9% | 21 |. That is, 
\C& g {mw ) — C' 8g (mw)\ can not satisfy the condition, larger than 2.4, because of the b —> s g constraint. Therefore, 
we come to the conclusion that the puzzle of B a —> can not get resolved within ler experimental bounds in the 

MSSM and model I and II 2HDMs. 

If one does not impose the b —> s g constraint it is possible to resolve the puzzle in the MSSM because the 
Wilson coefficient C 8g can reach values larger than 2.6 in some regions of parameter space. We have carried out 
detailed numerical calculations in the MSSM, imposing the constraints from the B ° — B% system, the mass difference 
A Md = (0.509 ± 0.004)ps _1 , mixing induced CP violation phase angle /? measured in charmonium B decays, sin 2/3 = 

however, without imposing the 


0.687 ± 0.032 Q, and B° —> JQy, in addition to the constraint from Bd 
b —> s g constraint. 8 F 8 ’ RR and 8 FR RL are constrained to be order of 10 -1 and 10” 2 respectively with moderate 
sparticle masses (say,500GeV) [23,123. In particular, Br (B° —> Xdj) < lx 10 -5 extracted from exclusive B —> p(w)y, 
as advocated in Ref. |2jj, gives a more stringent constraint. Br(R° —> Xd"f) directly constrains |C 7 7 (TOb)| 2 -|-|C 7 7 (TOb )| 2 
at the leading order. Due to the strong enhancement factor m g /mb associated with single S FRI ' RL ' 1 insertion term 
in Cy 7 (m&), 5 RR ( RL ' > (~ 10~ 2 ) are more severely constrained than 5 F 8 ^ RR \ However, if the left-right mixing of 
scalar bottom quark 5^ is large (~ 0.5), Sy^ RR> is constrained to be order of 10 2 since the double insertion term 
^ll(rr)^LR(LR*) a ] s0 eri i lallcec j by rn g /m b . In the case of LL (RR) insertion the precisely measured mass difference 
A Md imposes a more severe constraint on Sy 8 FRR . 

In numerical analysis we fix tan /3 = 10, vary m g and m g in the region between 300 GeV and 2 TeV, and the NHB 
masses in the ranges of 91GeV < mh ^ 135GeV, 91GeV ^ itih ^ 200GeV with m b < mn and 200GeV < to^ < 
240GeV for the fixed mixing angle 7t/2 of the CP even NHBs, and scan Sf 8 B in the range |< 5 i 3 lB | ^ 0.1 for A=B and 
0.05 for A ^ B (A = L, R). The numerical result for the correlation between branching ratios of B° —> 7r°7r° and 
B° K°K° is shown in Fig0 for the case of LL insertion. Due to the combined constraints mentioned above, in 
most of parameter space Cs g (mw) is not large enough to accommodate the data in ler region. However, there are 
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FIG. 3: The correlation between B° —> 7r°7r° and B° —> KK in MSSM with the LL insertion. 


small regions of parameter space with x 1 (say x ~ 40, x is the square of the ratio between scalar quark and gluino 
masses) where Cs g {inw ) is large enough to resolve the puzzle. In the case of both LL and RR insertion, the result is 
similar. In the cases of LR, both LR and RL insertions, we also have similar results. For x > 1 we have the regions 
with large enough Cg, g (mw) and the regions are larger than those in the cases of LL, both LL and RR insertions. In 
short, numerical results of Br show that the MSSM can explain the puzzle within ler experimental bounds under all 
relevant experimental constraints except that from Br (6 —> s g). 


Finally we consider the model III 2HDM (24| . In the model III 2HDM there are tree-level flavor changing neutral 
currents (FCNCL After diagonalizing the mass matrix of quark fields, the flavor changing (FC) part of the Yukawa 
Lagrangian is 1241 


= ^Qi,LH 2 U j>R + %Q ilL H 2 D jiR + h.c. 


D r 


(13) 


In order to obtain naturally small FCNC one assumes the Cheng-Sher parameterization 

rD x y/™?™] 

$ij ~ *ij v 

Phenomenological constraints on parameters of the models have been extensively discussed 


(14) 

For b —> dss, the 

couplings A bd,db,ss are involved. \m can reach 0.4 without a conflict with the measured mass difference AMs d if 


0 


the mass of pseudo-scalar Higgs boson Ma is large (say, ~ lTeV) m\ and it is also allowed by the recent data for 
B° —> p(u>)j. The constraint on A ss from the analysis on B° — B° s shows that A ss can reach 0(100) B which means 
that the coupling of Higgs to s quark is O(10 -2 ). It should be emphasized that the constraint from Bd —> P + P~ is 
irrelevant in the model III 2HDM because the decay involves A w besides A bd and A has no relation to A ss , which 
is different from the MSSM and 2HDMs I and II. 

In numerical calculations, we use = 120GeV, rrid = 6 MeV, Am = 0.3, \ ss = 150 and consequently obtain 
C 13 {mw) = 0.41. Corresponding to this value, \C% g (mw)\ ^ 0 . 6 , which can be satisfied under all relevant constraints. 
The numerical result for the values of parameters given above shows that 


Br(B u -> 7 r 7 r ) = 1.3 x 10 


-6 


Br(B {> -> A' u /\ u ) = 0.9 x 10 


That is, the data of branching ratios of B° —> Tr°n° and B° —> K°K° are accounted for, as expected. At the same 
time, we have checked that the NP contribution to Br(H° —> 7 r + 7 r“) is very small and negligible. 

In conclusion, we have analyzed the B° —> Tr°n° puzzle in three kinds of models beyond the SM. In the analysis lcr 
experimental bounds and 60% theoretical uncertainty which mainly comes from the input of non-perturbative param¬ 
eters have been taken into account. It is shown that the minimal flavor violation models, the minimal supersymmetric 
standard model, and the two Higgs doublet models I and II can not give an explanation of the B° —> 7 r° 7 r° puzzle 
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within lfj experimental bounds when all relevant experimental constraints are imposed and the model III 2HDM 
can explain the puzzle without a conflict with other experimental measurements. Therefore, if the data of Br for 
B° —> 7r°7r° and B° — > K°K° persist in the future, the MFV models, MSSM and model I, II 2HDMs will be excluded 
within lcr experimental bounds and the model III 2HDM will be survived to resolve the puzzle. As it is obvious, 
the analysis depends on the estimate of theoretical uncertainty. Our estimate comes from the uncertainties of input 
parameters (form factors, distribution amplitudes, CKM matrix elements, etc.) as well as the error estimate of ne¬ 
glecting electro-weak penguin and annihilation contributions. If the uncertainty were 70%, \Cs g (mw) — C' 8 g (mw)\min, 
which should be reached in a model in order to account for the data, would be 2.2 in the MFV models and 2.0 in 
the MSSM respectively so that the MFV models could not give an explanation of the data and the MSSM could. 
We have also analyzed the puzzle in the case without imposing the b —> s g constraint in the MSSM. For all kinds of 
insertions there are regions of parameter space where the puzzle can be resolved within lcr experimental bounds. We 
expect that similar effects appear in decays with PV final states, B —> 7r°p° and B —> K°K°*. 
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